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We conclude our analysis of the Nt = 6 equation of state (EOS) for two flavor QCD, first described at last 
year's conference. We have obtained new runs at am q — 0.025 and improved runs at am q = 0.0125. The results 
are extrapolated to m q = 0, and we extract the speed of sound as well. We also present evidence for a restoration 
of the SU(2)xSU(2) chiral symmetry just above the crossover, but not of the axial U(l) chiral symmetry. 



1. INTRODUCTION 

The equation of state (EOS) is essential for 
phenomcnological models of heavy-ion collision 
experiments that seek to detect the quark-gluon 
plasma (QGP). It is also of cosmological interest 
since the strongly interacting matter present in 
the early universe might have appeared in such a 
state. 

Following an initial EOS calculation at Nt = 4 
by some of us jl|, we reported preliminary re- 
sults for Nt — 6 in Melbourne Since then, 
we have completed runs at am q — 0.025 and ex- 
tended runs at am q = 0.0125 to smaller step size. 
With these additional simulations, we are able to 
fit our data to empirical expressions, based on 
standard sigma model assumptions with an 0(4) 
or mean field critical point at zero quark mass. 
The fits allow an extrapolation to am q = and 
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permit a smooth interpolation of the data, so, for 
example, the speed of sound can be determined. 
We have also extended the am q — 0.0125 data to 
higher temperature, roughly 250 MeV, where a 
plateau in the energy density is now evident. 

We refer the reader to Refs. |l],|| for details of 
our method. There we describe how the pressure 
and energy density are found from derivatives of 
the partition function with respect to the gauge 
coupling and bare quark mass and a nonpertur- 
bative beta function. 

2. SIMULATIONS 

The analysis requires both hot (12 3 x 6) and 
cold (12 4 ) lattices. Each hot (cold) simulation 
is at least 1800 (800) time units long after equi- 
libration. Near the crossover region the simula- 
tions were extended to more than 3000 units for 
the hot lattices. The range of gauge couplings 
and masses corresponds roughly to physical tern- 
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Figure 1. The EOS. The diamonds indicate an 
earlier result on N t = 4 lattices. 



peratures, < T < 250 MeV (based on the p 
mass 0), and mass ratios, 0.3 < m^/rrip < 0.7. 
We have also performed simulations at various R 
algorithm step sizes At to extrapolate results to 
At = 0. For instance, extrapolating the step size 
from 0.005-0.01 to zero produces a cumulative 
change in the pressure of 5% at 6/g 2 = 5.6. 

3. RESULTS 

In Fig. [j] we compare the Nt — 6 EOS with 
the Nt = 4 result and the continuum and lattice 
Stefan-Boltzmann laws. Here we have converted 
to a physical temperature scale. There is an ap- 
parent large finite size effect which is expected 
from the free lattice theory. For Nt = 6 we see 
that the results depend weakly on the quark mass. 
From the location of the maximum in the slope 
of (□) or ('tpip) with respect to 6/g 2 , we conclude 
that the pseudocritical temperature of the tran- 
sition is roughly 140 MeV (for both am q — 0.025 
and 0.0125). Figure [l| shows that the energy 
density is already substantial at this point. At 
high temperature the energy density has leveled 
off dramatically and its expected approach to the 
free lattice result is slow. 
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Figure 2. Fit to ipip an d its derivative with 
respect to 6/g 2 . The difference between 0(4) 
and mean field is not discernible with the present 
data. Extrapolation to m q — 0, however, gives a 
different critical coupling. 



To obtain a smooth interpolation of the above 
results and extrapolate to m q w (the physical 
value), we assume a second order phase transiton 
at wiq = (but see Ukawa's review ||] for caution- 
ary remarks). The nonanalytic part of the free 
energy obeys the scaling relation 

f crit (t,h) = b- d f CIit (byn,by»h) (i) 

where t and h, the usual scaling variables, are 
proportional to 6/g 2 — (6/<? 2 ) c and am q , respec- 
tively. Both the critical exponents, yh and yt, and 
the scaling function f cr it(t,h) are universal. Two 
flavor QCD has long been expected to be in the 
same universality class as the 3d 0(4) Heisenberg 
magnet [H. While the 0(4) critical exponents are 
known |]5[, the scaling function is not. Therefore 
we have performed new simulations of the 0(4) 
magnet to determine the scaling function. Re- 
cently it has been argued that the transition may 
be mean field || , for which the entire form of the 
scaling part of the free energy is known fij]. 

In Fig. U we show preliminary fits of ipip and 
its derivative to the mean field and 0(4) seal- 
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ing functions plus a polynomial in am q and 6/g 2 . 
For the quark masses used in our simulations, the 
fits are indistinguishable. However, the respective 
extrapolations to m q — are quite different. In 
Fig. H the mean field and 0(4) critical couplings 
correspond to zero mass temperatures T c ss 150 
and 160 MeV, respectively. When (□) and its 
derivative are included in the fits, these values are 
shifted down by about 10 MeV. The above indi- 
cates that present lattice simulations may still be 
too far from the scaling region and smaller quark 
masses are required to see the true scaling behav- 
ior. 

An extrapolation of the EOS to m q = is 
shown in Fig. || It is compared with the am q = 
0.0125 result, which reproduces the data reason- 
ably well. Note that here we fit (□), (M>tf>) and 
their derivatives with respect to 6/g 2 simultane- 
ously. The appearence of the bump in the energy 
density just after the transition is probably an 
artifact of the extrapolation (at m q = 0, the cor- 
responding region of 6/g 2 lies below the values 
of the coupling where we have done simulations). 
From Fig. || we again see a weak dependence on 
the quark mass. In Fig. |] we show the speed of 
sound calculated from this fit. The low temper- 
ature part of the curve is not trustworthy since 
the derivatives of the energy density and pressure 
are poorly known in this region. In fact, we ex- 
pect the hadron gas below the transition to have 
a nonzero speed of sound, which then dips down 
at the transition. 

4. CHIRAL SYMMETRY RESTORA- 
TION 

Sigma model analogies suggest that in the 
chiral limit of zero quark mass in the quark- 
gluon plasma phase, the SU(2)xSU(2) chiral 
symmetry is exact, but the anomalous U(1)a ax- 
ial symmetry may remain broken Q). Early ef- 
forts to assess the status of these symmetries 
looked for the expected degenerate multiplets in 
screening masses. However, because these stud- 
ies did not include difficult-to-measure quark-line 
disconnected graphs, they were inconclusive ||. 
The large data sample generated in the EOS 
study makes it feasible to reexamine this ques- 



Fit using 0(4) scaling 




0.10 0.15 0.20 0.25 

T(GeV) 



Figure 3. The EOS extrapolated to am q = 0.0. 
The fit and data at am q = 0.0125 are also shown. 
The triplet of fit curves at both masses, as la- 
beled, indicates the central value and one stan- 
dard deviation spread resulting from the statisti- 
cal uncertainty in the fit. 

tion. Knowing the answer gives important insight 
into the mechanics of the phase transition. The 
Columbia group has also pursued this question 
0, the Bielefeld group has reported preliminary 
results |l0| , as has Lagae [0 

To address this question we focused on suscep- 
tibilities for the following hadron channels: /o 
(also known as a or iptp), o,q (5), and 7r. For exam- 
ple the /o susceptibility is constructed from the 
point-source correlator as follows: 

X/o = J d 4 x(/ (x)/o(0)) - (/o(0)) 2 (2) 

From these susceptibilities, we define the order 
parameters 

XU(1) = Xtt - Xao', XSU(2) = X/o - Xtv- (3) 

If the order parameter is nonzero in the chiral 
limit, the corresponding symmetry is broken. 

The /o susceptibility can be decomposed into 
quark-line-connected and disconnected parts as 
X/o = Xconn + Xdisc- The connected part is iden- 
tified with Xa 1 so is obtained from the point 
hadron correlator using traditional methods. The 
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Figure 4. The speed of sound for am q = 0.0125. 
The three curves indicate the central value and 
the one standard deviation spread resulting from 
the statistical uncertainty in the fit. 



disconnected part is the variance from configu- 
ration to configuration of the space-time volume 
average of /o = ipip. The latter is conventionally 
measured from Tr(p+m g ) _1 , which we estimated 
using multiple random sources on each configura- 
tion. Finally, the pion susceptibility Xic can a l so 
be obtained from (fa)/m. 

For this study we fixed 6/g 2 = 5.45 and varied 
the quark mass from am q = 0.0075 to 0.025, cor- 
responding to the high temperature phase. The 
crossover to the low temperature phase occurs 
at a slightly higher mass, and the chiral limit 
corresponds to a temperature of approximately 
1.2T c (ma = 0). 

Preliminary results are shown in Fig. ^. In the 
two-flavor staggered fermion scheme the fermion 
determinant is rigorously even in the quark mass. 
In consequence the order parameters are also even 
in am q . For T > T c there are no known infrared 
singularities, even in the chiral limit. Thus we 
have indicated a linear extrapolation in (am,) 2 , 
leading to a result consistent with the sigma 
model scenario: a restoration of SU(2) x SU{2) 
but not of U{1)a (approximately 2a). 

This work was supported by the US DOE and 
NSF. Computations were done at the San Diego 
Supercomputer Center, the Cornell Theory Cen- 
ter, Indiana University, and the University of 
Utah. 
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